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The diffusive motion of a colloidal particle trapped inside a small cavity filled with fluid is reduced 
by hydrodynamic interactions with the confining walls. In this work, we study these wall effects on 
a spherical particle entrapped in a closed cylinder. We calculate the diffusion coefficient along the 
radial, azimuthal and axial direction for different particle positions. At all locations the diffusion 
is smaller than in a bulk fluid and it becomes anisotropic near the container's walls. We present a 
simple model which reasonably well desribes the simulation results for the given dimensions of the 
cylinder, which are taken from recent experimental work. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

Modeling the motion of colloidal particles under con- 
finement is of interest for many industrial applications 
which involve the entrapment of particles in microporous 
materials and diffusion controlled reaction mechanisms. 
An important class is that of coating processes where the 
particles have to diffuse towards a wall before being im- 
mobilized. Other applications involve microfiuidic chips 
P, [|[ and surface-based biosensor Q . Furthermore the 
confinement of particles in more than one dimension is 
relevant for the synthesis of colloids inside droplets [4[ 
and for understanding the transport of proteins, lipid 
granules and other macromolecular assemblies through a 
biological cell. 

Most of the past work on modeling the mobility of a 
particle in confinement has focussed on a particle in prox- 
imity of a flat wall or within an infinitely long cylinder 
|5H9j . The motion of an arbitrarily shaped particle near 
boundaries has been studied in , while the settling of 
a particle inside a cylinder closed at the bottom but char- 
acterized by a free surface at the top has been analyzed 
in (Tlj . Finally for a historical perspective on the mo- 
tion of a sphere in an uncompressible fluid, from Stokes' 
law to Faxen's solution of the Navier-Stokes equations 
for particular bounded systems, see also [12|. Under the 
assumption of creeping motion, previous investigations 
are mainly based on various approximate perturbativc 
schemes, which are generally not valid when the size of 
the particle, the size of the cavity and the gap between 
particle and wall are all comparable, which is the case we 
are interested in. Recently, however, a semi-analytical 
procedure for the calculation of the mobility matrix of a 
spherical particle inside a cylinder has been presented in 
13], which can be used under a wide range of confinc- 
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ment conditions. 

The motion of a particle in a closed three dimensional 
container has been studied less intensely. The flow field 
around a sphere has been studied in [3 [HI m vari- 
ous confined geometries. The motion of a particle along 
the axis of closed cylinders and closed cones is discussed 
in ref. [16j |. where experimental and theoretical results 
are presented for a millimeter sized sphere settling in a 
very viscous silicone oil. Agreement is found until the 
wall-particle separation is of the order of the particle di- 
ameter. Recently the diffusion coefficients of /im sized 
spheres moving in a squat cylinder containing a water 
solution has been studied [l?} , showing that the diffusion 
is strongly affected by the confinement and substantial 
corrections to the Stokes' law are required. 

The analysis of the motion of a particle inside a cavity 
can be based upon the solution of Navier-Stokes equa- 
tions under particular regimes, such as the creeping mo- 
tion. However, such equations are not always easily 
solved, especially for complex geometries or when there 
are many particles in suspension. Moreover they are 
based on the hypothesis that the solvent is a continuum. 
But in conditions relevant for microfiuidic devices, this 
hypothesis is not always fulfilled and fluctuations of the 
solvent velocities may be relevant. To study diffusion 
of colloids in microcavities it is therefore important to 
test different tools, like simulations in which the solvent 
is explicitly introduced, or in which the solvent velocity 
is allowed to fluctuate. In this work we use the Multi- 
Particle Collision Dynamics technique, initially proposed 
by Malevanets and Kapral (is) . It consists of a coarse- 
grained model to describe the solvent completed with 
specific rules (see below) to describe solvent-solute inter- 
actions. It has been successfully used to study hydrody- 
namic interactions in complex fluids, at equilibrium and 
far from it. For a review, see (l9j . 

We will briefly introduce the simulation technique in 
section |TTJ In section Hill we discuss the measurements of 
friction and diffusivity. In section HVl we will analyze the 
simulations results for the diffusion of a spherical particle 



2 



inside a microcavity in the shape of a closed cylinder. We 
will show that the friction is different along the radial, 
azimuthal and axial directions, and that the correspond- 
ing diffusion coefficients change in different ways with the 
particle's radial and axial position. In section [V] we will 
compare our results with experimental and theoretical 
data. A simplified model to interpret our results is given 
in section IVTl Finally, we give our conclusions in section 

Em 



II. SIMULATION TECHNIQUE 

A severe difficulty in simulating colloidal solutions is 
to correctly deal with the large range of time and length 
scales involved. Whereas the solvent velocities generally 
relax in less than a picosecond, colloids typically diffuse 
over a length equal to their own size in seconds. For- 
tunately it is often not necessary to exactly reproduce 
such time and length scales. A coarse-grained approach 
is possible in which the scales are telescoped together [2(| , 
while keeping an order of magnitude difference between 
each relevant time scale. 

As already mentioned, we will use Multi-Particle Col- 
lision Dynamics (MPCD) pj, ^ to describe the solvent 
degrees of freedom, with no-slip boundary conditions for 
the solvent-solute interactions. The implementation of 
the no-slip boundary conditions has been discussed in 
detail in |21j . The dynamics of the system is made up of 
two steps: streaming and collision. In the streaming step, 
the position and velocity of each particle is propagated 
for a time St by solving Newton's equations of motion. 
During the collision a lattice is placed in the system [22j . 
This lattice is typically cubic, with a lattice constant a. 
Solvent particles are then attributed to the cells in which 
they happen to be positioned. Next the velocities Vi of all 
particles are stochastically rotated, relative to the center 
of mass of the corresponding cells according to the for- 
mula: 

v' i (t)=u + Sl{v i (t)-u}, (1) 

where u is the mean velocity of the particles within the 
collision cell to which particle i belongs and l~i is a ma- 
trix which rotates velocities by a fixed angle a around 
a randomly oriented axis. Through the stochastic rota- 
tion of the velocities, the solvent particles can effictively 
exchange momentum without the need of introducing di- 
rect forces between them. As the collision step preserves 
mass, linear momentum and energy, the correct hydro- 
dynamic behavior is obtained on the mesoscopic scale 

[31. 

When colloids are present, Newton's equations of mo- 
tion are solved also for them during the streaming step. 
However special care must be taken when no-slip bound- 
ary conditions are applied on the colloidal surface and 
on the confining walls. The implementation of the no- 
slip boundary conditions is described in the following. 



Streaming step: when a solvent particle crosses the col- 
loid (or wall) surface, it is moved back to the impact 
position. Then a new velocity is extracted from the fol- 
lowing distributions for the tangent (v t ) and the normal 
component (v n ) of the velocity, with respect to the sur- 
face velocity: 

^-^T eXP ("2l^J' (2) 

, . / m / mv? \ 

Here m is the mass of the solvent particle, ks the Boltz- 
mann's constant, and T the temperature of the system. 
Once the velocity has been updated, the particle is dis- 
placed for the remaining part of the integration time step. 
Collision step: virtual particles (VP) are inserted in those 
parts of the collision cells which are physically occupied 
by the colloid or by the container walls. The VP density 
matches that of the MPCD solvent, while their velocities 
v( p are obtained from a Maxwell-Boltzmann distribu- 
tion, whose average velocity is equal to the local velocity 
of the colloid or to the velocity of the walls, and the tem- 
perature is the same as in the solvent. Notice that we 
do not take rotations of the particles into account. VP 
are sorted into the collision cells according to their coor- 
dinates, exactly as the MPCD particles. This method is 
similar to that used in (23 - [27j . Then the average velocity 
of the center of mass of the cell is computed as: 

umpcd nvp 

u _ E^w+E^r^ (4) 

n-MPCD + nyp 

where umpcd and nyp are the number of MPCD par- 
ticles and VP, respectively, belonging to the same cell. 
Velocities of both MPCD and VP belonging to the same 
cell are rotated according to the rule given in eq.Q. 

Due to the exchange of momentum between the sol- 
vent and the colloidal particle, the force exerted upon 
the latter can be expressed as F = f s + f c , where f s and 
f c are the forces during the streaming and the collision 
step, respectively. The former can be calculated as: 

i— 1 

where Q is the number of MPCD particles which have 
crossed the surface of the colloid between two collision 
steps, and Ap^ = mAvi is the change in momentum of 
solvent particle i, which has been scattered by the colloid. 
The force exerted during the collision step is: 

i— 1 

where q is the total number of virtual particles which 
belong to a tagged colloid, and Ap^ p is the change in 
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momentum of the virtual particle i as a consequence of 
the collision step. 

The force exerted by the solvent on the colloidal par- 
ticle, will be used in section IIVI in order to measure the 
friction on the particle itself. 

In our simulations we use a collision angle a = n/2, 
an average solvent particle density of 5 particles per col- 
lision cell, and a collision time step 5t = O.lto, where 
to = ay/ m/ksT is the simulation unit of time. This 
choice leads to a mean-free path which is an order of 
magnitude smaller than the collision cell size a, which is 
important to simulate the dynamics of a colloid in a liquid 
instead of a gas [20[. Simulations have been performed 
with a colloidal sphere of radius R co i = 3a, which is suffi- 
ciently large to accurately resolve the hydrodynamic field 
to distances as small as i? co //6, as already shown in Refs. 
[2M [H| . The radius and height of the closed cylinder are 
Rcyi = H cy i = 54a. In this way we can make a direct 
comparison with the experimental case discussed in [l7| . 
The typical length of a simulation run is about 1.5 • 10 6 
integration time steps. 



III. MEASURING FRICTION AND 
DIFFUSIVITY 

The drag force F and torque T acting on a particle 
that moves through a fluid with a steady velocity V and 
steady angular velocity £1 can be expressed as 



M 



(7) 



where M is the friction tensor represented by a 6 x 6 
matrix. The friction matrix can be subdived in 4 blocks, 



7* ctr 
M = ( " t " 



(8) 



where the superscript t stands for translation, r for rota- 
tion and tr and rt for the coupling between rotational mo- 
tion and (translational) force and between translational 
motion and (rotational) torque, respectively. The blocks 
of the matrix M can be obtained in the following man- 
ner. During the simulation the colloidal sphere is kept in 
a fixed position and fixed orientation by applying an ex- 
ternal constraint force F c and external constraint torque 
T c . The constraint force and torque are equal, but op- 
posite, to the total force and torque exerted by the sol- 
vent on the sphere F c = — F and T c = T. Through 
the fluctuation-dissipation theorem the autocorrelation 
functions of F c and T° can be connected to each of the 
components of the friction tensor. For example the ex- 
plicit expression for S* [H, [2t| is: 



77* 

"a/3 



1 



ksT jo 

Similarly we can also measure 
1 



"a/3 



dr(F Q c (r + r)F|(r )) ro . (9) 
dr(F^T +r)TI(r )) T . (10) 



FIG. 1: For a generic sphere position, we plot the unit vectors of 
the radial (r), azimuthal (0) and axial (2) directions. The distance 
from the cylindrical walls is <5, while the distance from the bottom 
wall is z. When the sphere is in the mid-plane z = H cy i / 2, when 
the sphere is placed along the axis 8 = R cy i- 



and so on. The quantities a, j3 refer to two of the cylin- 
drical coordinates r, 9, z and the average is taken over 
many different time origins To. As shown in Fig. (JT|), for 
a given position of the sphere we can define an orthog- 
onal frame with unit vectors pointing in the local r, 9 
and z directions. Because the sphere is in a fixed posi- 
tion, this orthogonal frame does not change during the 
simulations. 

The friction tensor M is connected to the diffusion 
tensor through the Einstein relation D = ksTM~ 1 . The 
diffusion tensor can be subdived in blocks too: 



D 



(11) 



where D t is the translational diffusion tensor, D r the ro- 
tational diffusion tensor, and D tr and D rt determine the 
amount of cross-correlation between translational and ro- 
tational displacements. 

A strong coupling between translation and rotation is 
expected when the particle gets very close to a wall. This 
will cause a torque- free particle to rotate when it is trans- 
lated along the wall, and vice versa a force-free particle 
to translate when it is rotated near a wall. However, 
after a preliminary analysis on the walls' effects, based 
upon a sphere nearby a flat wall Q and inside an in- 
finitely long cylinder Q, for the particular gaps we are 
interested in, we find that such a coupling is very small: 
in the worst case scenario (d/R — 0.5) the effect of in- 
cluding translation-rotation coupling on the translational 
diffusion is less than 2%. We expect that the contribution 
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FIG. 2: Running integral of the autocorrelation function of the 
constraint force F c applied to the sphere in the mid-plane. The 
distance S from the cylindrical wall is in units of the particle diam- 
eter. Time is in units of to = a{m/ kgT) 1 / 2 and friction in units 
of m/trj. Upper panel: radial friction. Bottom panel: azimuthal 
friction, which is very similar to the axial friction (not shown in 
this picture). 



due to the coupling between rotation and translation is 
small in the closed cylinder as well, so that we can neglect 
the terms 5 tr and S rt in the friction tensor. 

In this paper we focus on two cases: one in which the 
particle is set in the mid-plane of the cylinder and an- 
other in which it is located along its axis. For both these 
cases, symmetry requires that 3* is diagonal. Therefore 
the translational diffusion tensor D* is also diagonal and, 
under the hypothesis that the coupling between rotation 
and translation is negligible, we set D* rr = fcBT/3* r , 
D gg = ksT/Sgg and D\ = k B T/~} zz . 

Because in this paper we measure only the transla- 
tional diffusion and not the rotational diffusion, we drop 
the superscript t and simplify the notation as D r , Dg and 
D z . Some specific components of the running integral of 
the autocorrelation function of F c are shown in Fig. ©: 



FIG. 3: Diffusion coefficients along the radial (r), azimuthal (8) 
and axial (z) direction. Upper panel: the sphere is in the mid- 
plane and 8 is the separation from the cylindrical walls. Bottom 
panel: the sphere is placed along the cylinder axis and z is the 
distance from the closest end- wall. Distances are in units of the 
particle diameter. Diffusion is in units of the bulk value Do (see 
main text). 



the sphere is in the mid-plane but at different distances 
S from the cylindrical walls. In particular in the upper 
panel we show the radial component of the friction 5 rr , 
while in the bottom panel the azimuthal one Sgg. 

Finally we note a peculiarity which is linked to the 
fact that in our simulations a colloidal particle collides 
with only a few solvent particles within each particular 
time step St. If we consider a particular direction a, 
there are essentially two contributions to the total fric- 
tion: one coming from the local Brownian collisions with 
the particles of the fluid (£ B ), while the other is due 
to the long-range hydrodynamic interactions (flow fields) 
(£ H )- A simple empirical formula says that the hydro- 
dynamic and the Brownian friction should be added in 
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parallel j3(J HH in order to obtain the total friction: 

i/E aa = i/e + i/ef Q (12) 

£ B does not depend on the sphere's position and corre- 
sponds to the height of the short-time peak of the integral 
in eq. ([9]); its value is 645 ± 1 (in units m/to) for all di- 
agonal components. The long-time limit of the profiles 
shown in Fig. ([2]) provides the total friction upon the 
sphere. Thus, the hydrodynamic term can be extracted 
by inverting ea. (IT2l . For validation of this method and 
further applications of this concept see [2(| Hll ■ 

IV. SIMULATION RESULTS 

The measured diffusion coefficients are summarized in 
Fig. Q for the mid-plane (upper panel) and for the cylin- 
der axis (bottom panel). The diffusion is normalized to 
the bulk value Dq for a sphere with no-slip boundary 
conditions upon its surface. The latter is provided by 
the Stokes-Einstein law: Dq = ksT / (Q j nvR co i), where v 
is the solvent viscosity. 

First we focus on the case where the sphere is located in 
the mid-plane of the cavity. When the colloidal sphere is 
set near the center of this mid-plane (large 5), we expect 
that D r = Dg because of symmetry. Actually we ob- 
serve that all three diagonal components are almost the 
same and their value is significantly lower (25%) than in 
a bulk system. When the sphere is shifted towards the 
cylindrical wall, the diffusion decreases even more. For 
off-centered positions, the radial component D r decreases 
more rapidly than Dg and D z do. 

Next we focus on the case where the sphere is located 
on the central axis of the cavity, at a distance z from the 
nearest end-wall. Although in this case technically no 
azimuthal direction can be defined, we have made two 
orthogonal measurements perpendicular to the z-axis, 
and labeled them as D r and Dg in the bottom panel 
of Fig. ([3]). As expected, because of symmetry we find 
D r = Dg within uncertainties. Close to the end-wall, D z 
decreases much faster than the other two components do. 
This situation is similar to the case of a sphere close to a 
flat wall. However the similarity is only qualitative and 
not quantitative, as we will discuss in section fVTl 

V. COMPARISON WITH EXPERIMENTAL 
AND THEORETICAL PREDICTIONS FOR THE 
CLOSED CYLINDER 

First we compare our data with some theoretical re- 
sults published recently and next with some recent exper- 
imental results. It is interesting to notice that our results 
presented in the upper panel of Fig. © are in agreement 
with the findings in et al. (l3j . where hydrodynamic cal- 
culations are performed of the drag on a particle in an 
infinitely long cylinder. In particular the authors observe 
that for R cy i > 4i? co ;, all three diagonal components of 



TABLE I: Axial friction upon a sphere in a closed cylinder. The 
radius and the height of the cylinder are in units of the sphere 
radius. Theoretical predictions corresponds to Il6l . The friction is 
normalized to the bulk value which is £o = 67rfR co ;. 
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the friction increase monotonically when the particle ap- 
proaches the wall. The radial friction is expected to vary 
stronger than the other two frictions, because its motion 
is perpendicular to the confining walls, while in the other 
cases the motion is tangent to the wall. As discussed in 
[l3j , the resistance for the radial motion is inversely pro- 
portional to the gap between particle and wall, while for 
the tangential motions it varies logarithmically with the 
gap. Qualitatively our results are similar, in the sense 
that the azimuthal diffusion for small particle-wall sepa- 
rations varies slower than the radial diffusion does; this 
suggests that the role of the cylindrical wall is dominant 
in this region. However, we will show that for our squat 
cylinder the role of the end walls cannot be neglected at 
any sphere position (section EH- 

We have performed two additional simulations in order 
to compare the results with the theoretical predictions in 
16] for the axial friction j z = £z/£o in a closed cylin- 
der. Here £ = 6nvR co i corresponds to the friction on a 
sphere in an unbounded system. The two cases that we 
have chosen are, one in which the cylinder is very narrow 
and long, and one in which it is squat. The results are 
summarized in Table |U The agreement is good. 

In case the closed cylinder is very long, the friction on a 
sphere in the mid-plane should not be very different from 
that predicted for a sphere in an infinitely long tube. For 
example, for the case R cy i — ^R co i the friction inside an 
infinitely long cylinder can be found in [13[ . In particular 
for a sphere placed on the cylinder axis, the values are: 
j r = 1.763 and j z = 1.992. Our simulation results for 
a closed cylinder with the same radius but finite length 
H cy i = 20R co i agree within error bars: j r — 1.79 ± 0.07 
and -f z = 2.03 ± 0.08. We can conclude that, for a long 
cylinder the role of the end walls is minimal. 

The dependence of the diffusion on the spatial posi- 
tion of a particle inside a closed cylinder has been re- 
cently studied through experiments in [l7|. The authors 
discuss the case of a \im size particle in an aqueous sus- 
pension. The microcavity is produced via soft lithogra- 
phy where PDMS membranes are used. In Fig. ((4]) we 
compare their experimental results (open symbols) with 
our simulation results (filled symbols) in the mid-plane 
for the radial and the azimuthal components. The setup 
of the experiments did not allow the measurement of the 
axial diffusion. The agreement between the two sets of 
data is good. In particular the region where D r and Dg 
begin to differ is well captured. 
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FIG. 4: A comparison between simulation results (filled sym- 
bols) and experimental data [17j (open symbols) for the diffusion 
coefficient in the mid plane of a closed cylinder. Circles: radial dif- 
fusion. Squares: azimuthal diffusion. The wall-particle distance 8 
is in units of the particle diameter d. R cy i = 9d and H cy i = 9d. For 
sake of clarity, we show the error bars for only a few experimental 
data points. 



The good agreement with both experimental and the- 
oretical predictions reaffirms that MPCD simulations en- 
able us to predict quantitatively the friction on a colloidal 
particle in a confined space. 



VI. DISCUSSION OF WALL EFFECTS 

We will now discuss the effects of confinement on the 
diffusivity of a colloidal particle, by focussing on the ef- 
fects of each individual type of boundary, i.e. cylinder or 
flat wall. 

Generally speaking, for small Reynolds numbers, the 
correction term to the friction on a confined particle is 
7 = £/£o: where £o is the friction on a sphere in an 
unbounded solvent. In principle 7 depend on the direc- 
tion of motion. Since there are no sufficiently general 
solutions available of the Navier- Stokes equations for a 
closed cylinder, we will try to estimate the total wall cor- 
rective factor as a combination of the corrective terms 
due to simplified geometries, such as the infinitely long 
cylinder, flat walls and square ducts. For ease of reference 
we reproduce explicit formulas available in the literature. 

The wall correction for the axial friction on a sphere 
at an eccentric position is, approximately, [5[: 



I? 1 = 1 



2.10444-0.6977 




This expression is valid mainly near the axis of the cylin- 
der. If the sphere is exactly on the axis, a more precise 



FIG. 5: Axial diffusion coefficient in the mid-plane. Simula- 
tions (diamond symbols) are compared with: i) prediction for an 
infinitely long cylinder, valid near the axis (dashed line), ii) pre- 
dictions for a flat wall, valid at the periphery of the cylinder (solid 
line), and iii) predictions from numerical solutions of the Navier- 
Stokes equations in an infinitely long cylinder (dotted line). All 
these predictions disregard the role of the end-walls, so that the 
diffusion is globally overestimated relative to our simulation re- 
sults. 



correction due to Bohlin [32j, also reported in Q, is: 



v c V l 



[1 - 2.10444x+ 2.08877x 3 - 6.94813x 



5_ 
-1 



-I.372x 6 + 3.87x 8 - 4.19x 10 ] , (14) 

where \ = R co i/Rcyi- When the sphere is small (relative 
to the cylinder's diameter) and very close to the cylinder 
wall, it is expected that the curvature of the wall is not 
very important and that the sphere behaves like near to 
a flat wall [33| ■ In such a case the wall correction for the 
parallel motion is @: 



R 



1 16 V h 
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(15) 



while for the perpendicular motion it is: 



7_l 



al1 = -sinh(a)£. 



n(n + 1 



3 ' K ~' ^ (2n-l)(2n + 3) 

2 sinh(2n + l)a + (2n + 1) sinh 2a 
4sinh 2 (rc. + 0.5)a - (2n + l) 2 sinh 2 (a) 



- 1 



(16) 

Here a = cosh -1 (h/R co i), in which h is the distance of 
the sphere from the wall. 

In Fig. ([S]) we plot the simulation data for D z and 
the theoretical predictions for an infinitely long cylin- 
der. We use eq. (|13[) which is valid near the cylinder axis 
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FIG. 6: Rectangular parallelepiped circumscribing the cylinder. 
For a generic particle position, the indexes 1 and 1 ' refer to the walls 
perpendicular to the radial motion; 2 and 2' to the walls perpendic- 
ular to the azimuthal motion; 3 and 3' for the walls perpendicular 
to the axial motion. 

(dashed line) and eq. (IT5t which is valid near the the 
cylinder walls (solid line). The profile provided for the 
region nearby the cylindrical wall is also in good agree- 
ment with the model in [8]. Although such equations 
have only a limited range of validity, in Fig. ([5]) they have 
been plotted along the entire range of distances we are 
interested in. In the same picture we also show the diffu- 
sion obtained from the numerical solution of the Navier- 
Stokes equations (dotted line) for a cylinder provided by 
S. Bhattacharya (private communication), according to 
the method described in [13J. We see that the latter 
result agrees very well with eqs (fTBI and (|T5j) in the ap- 
propriate limits. The discrepancy with our simulation 
results, however, is quite large for all values of S. This 
comparison clearly shows that the end walls' effects are 
important for all particle positions in the mid-plane of 
the cylinder both close to the axis, as well as in proxim- 
ity of the cylinder walls, in case the separation between 
the end- walls is of the order of the cylinder diameter or 
less. 

In the following we will attempt to find a simple model 
that can predict the combined effects of all walls, at least 
with the dimensions of our simulation box. 

It has been shown [34| that when the particle is placed 
on the axis of a container and the size of the particle is 
small compared to the width of the container, the ax- 
ial friction depends only weakly on the precise geometry, 
whether cylindrical or rectangular. Inspired by this ob- 
servation, we approximate the closed cylinder by a rect- 
angular parallelepiped circumscribed to it, as shown in 
Fig. ©. We then express the total correction factor T 
for the friction as a combination of single flat-wall cor- 
rections 7: 

T = 7i( rf i) + 7v{di>) + 72(^2) + 72' (d 2 <) 

+73(d3)+73'(«*3')-(n-l), (17) 

where di refers to the distance of the particle from flat 
wall "i", as numbered in Fig. ((5]), n is the number of 



boundaries or walls included in the model. In the par- 
allelepiped model, n — 6. We have to subtract n — 1 
in order to avoid overcounting the bulk value; this is a 
mere consequence of the definition of the wall correction 
factor. 

Before continuing with the closed cylinder, let us con- 
sider how our model performs for infinitely long systems, 
in particular for an infinitely long cylinder and for an 
infinitely long square duct. Our model in this case has 
just for flat walls, since we must remove the walls 3 and 
3' in Fig. ©; consequently n = 4 in this case. For the 
cylinder we have j ax is — 1.132, when R cy i = 18i? co z; for 
the square duct we have 7^ = 1.118 when the half- width 
of the square is equal to 18i? co ; [34[ and for our model 
we find r = 1.129. All three numbers agree reasonably 
well. 

Coming back to the case of the closed cavity, the mo- 
tion along the radial direction can be decomposed into a 
motion perpendicular to two walls and parallel with re- 
spect to the remaining four walls. Similar conclusions can 
be drawn for the azimuthal and axial motion. So we can 
specialize the relation in (|17l) in the following manner: 

r r = 7±,l( rf l) +7-L,l'(dl') +7||,2(^2) +7||,2'(^2') 

+711,3(^3) +7||,3'(d3')- 5 (18) 
r = 7||,i( d i) +7||,i'(di') +7^,2(^2) +7-L,2'(e?2') 

+7||,3(«fe) +7||,3'( rf 3')- 5 (19) 
T z = 7||, i(di) +7||,i'(^i') +711,2(^2) + 7||,2'(d2') 

+7x,3(d 3 )+7J.,3'(d8')-5 (20) 

where the expression for the parallel (||) and perpendic- 
ular (_L) corrections are given by eqs (fl~5|) and (fTB]) . re- 
spectively. In Fig. |(7J) we compare these predictions with 
our simulation results; the upper panel refers to the mid- 
plane, while the bottom panel to the particle placed along 
the container axis. The agreement is very good, which is 
rather surprising given the simplicity of the model. Note 
however, that the radial and azimuthal diffusion coeffi- 
cients predicted by our model are not exactly equal on 
the axis, as they should be. Its success is possibly caused 
by a subtle cancellation of errors, which we will discuss 
next. 

The rectangular parallelepiped model is based upon a 
linear superposition of the walls' effects. In this approx- 
imation, we disregard the influence of the curvature of 
the vertical walls and that of the corners obtained by the 
intersection of the cylindrical wall and the top and bot- 
tom walls. For a particle in the mid-plane and close to 
the vertical wall, ignoring the effect of the curvature, the 
presence of the corners obtained by the intersection of 
the cylindrical wall and the top and bottom walls leads 
to larger frictions, so that their neglect leads to an under- 
estimation of the friction and consequently to an overes- 
timation of the diffusion coefficients. 

On the other hand, using the superposition model for 
the top and bottom walls per se leads to a slight overes- 
timation of the perpendicular frictions and thus to slight 
underestimations of the corresponding diffusions. For the 
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FIG. 7: Upper panel: diffusion in the mid-plane, where <5 is the 
distance from the cylinder wall. Bottom panel: diffusion along the 
axis, where z is the distance from the bottom wall. Symbols are 
for simulation data. Lines for the rectangular parallelepiped model: 
D r (dashed line), Dg (solid line), D z (dotted line). 



slit geometry this may be seen in Fig. (4) in ref.[35|. The 
magnitude of such effects depends on the distance of the 
particle from the walls. For the present case, the parallel 
friction for a slit geometry is quite well represented by 
the linear superposition model, whereas the perpendicu- 
lar friction is overestimated. 

For the configurations studied here the overestimation 
of the friction perpendicular to the end walls is large 
enough to compensate the underestimation of the total 
friction which is intrinsic to the superposition model ap- 
plied to a geometry with corners. 

The above discussion can also be cast in the frame- 
work of the method of reflections, where the role of walls 
is mimicked through a set of images of the colloidal par- 
ticle. This set of images is used to enforce a zero velocity 
boundary condition for the fluid near the planar surfaces. 



FIG. 8: Comparison between simulation data and the superposi- 
tion model applied to a combination of an infinitely long cylinder 
and two fiat walls. Symbols are for simulation data. Lines for the 
rectangular parallelepiped model: D r (dashed line), Dg (solid line), 
D z (dotted line). 



For a real parallelpid the superposition model only takes 
into account the six images due to the single reflections 
on each of the sides of the parallelepiped. However, ev- 
ery time a new image is introduced (to treat a new wall) 
the fluid velocity profile on the existing walls is modi- 
fied. This is then corrected, usually in an iterative fash- 
ion. Corner effects arise from the reflection of each image 
particle with respect to two perpendicular walls. In the 
slit geometry, instead, successive corrections to the linear 
superposition, comes from the multiple reflections of the 
images, similarly to those that appear in a barber shop. 

We have also tried a superposition of the friction ex- 
pressions for an infinitely long cylinder (private commu- 
nication from Bhattacharya) and two planar walls. In 
this case the agreement with the simulation data is less 
good, especially for the radial and azimuthal diffusion as 
shown in figE] This result shows that the superposition 
model is quite sensitive to the details of the surfaces we 
apply it, and therefore it should be applied with care. 



VII. CONCLUSIONS 

We have studied, by means of particle based simula- 
tions, the diffusion of a spherical particle inside a micro- 
cavity in the shape of a squat closed cylinder. We have 
analyzed the diffusion in the middle plane as a function 
of the particle-cylinder separation S and along the axis 
of the cylinder as a function of the particle-end wall sep- 
aration z. For our squat cylinder we find that when the 
sphere is near the center of the cylinder the diffusion is 
the same along the radial, azimuthal and axial direction 
and it is significantly lower than in the bulk. In the mid- 
plane, when the particle is shifted towards the cylinder 
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wall, the diffusion decreases more in the radial direction 
than in the azimuthal and axial ones. The last two are 
very similar to each other. When the particle is shifted 
along the axis we find that the diffusion decreases much 
more along the axial direction than perpendicular to it. 
In general we can say that the diffusion in the direction 
mostly parallel to the closest confining wall, decreases 
slower than in the direction perpendicular to it. 

We have compared our simulations with experimental 
results for a /im sized particle in a cylindrical microcav- 
ity, finding good agreement. We also find good agreement 
with the theoretical predictions for a closed cylinder re- 
ported in [16j. 

Finally we have shown that the simplest approximation 
to capture confinement effects, namely one which uses the 
result for an infinitely long cylinder for the axial diffusion, 
does not work. Neither does treating the closest wall as 
a flat wall when the particle is at the periphery of the 



box. The end-wall effects appear very important for all 
positions of the sphere inside the squat cylinder. We 
have investigated a simplified model in which the effect 
of all walls is approximated as a superposition of the 
effects of individual walls of a rectangular parallelepiped 
circumscribing the cylinder. This model works well for 
diffusion coefficients both in the mid-plane and along the 
axis, but it owes its success to a fortunate cancellation of 
errors. 
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